Fluctuations relation and external thermostats: 
an application to granular materials 
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In this note we discuss a paradigmatic example of interacting particles subject to non conservative 
external forces and to the action of thermostats consisting of external (finite) reservoirs of particles. 
We then consider a model of granular materials of interest for experimental tests that had recently 
attracted lot of attentions. This model can be reduced to the previously discussed example under 
a number of assumptions, in particular that inelasticity due to internal collisions can be neglected 
for the purpose of measuring the large deviation functional for entropy production rate. We show 
that if the restitution coefficient in the granular material model is close to one, then the required 
assuptions are verified on a specific time scale and we predict a fluctuation relation for the entropy 
production rate measured on the same time scale. 
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I. NONEQUILIBRIUM 

In studying stationary states in nonequilibrium sta- 
tistical mechanics, [1, 2], it is common to consider d- 
dimensional systems of particles in a (finite) container So 
forced by non conservative forces whose work is controlled 
by thermostats consisting of particles moving outside So 
and interacting with the particles of So through interac- 
tions across the walls of S , [3]. If X° = (x?, . . ., Xjy ) 
are the particles positions in an inertial system of coor- 
dinates, the equations of motion are determined by their 
mass m, by the potential energy of interaction V^X ), by 
external nonconservative forces Fi(X ,^) and by ther- 
mostat forces — i?j as 

m x ) ° = -3 x .y(X°)+F i (X () ;$)- ) ? i , (1) 

where i = 1, . . . , Nq and # = (tpi, . . . , ip q ) are strength 
parameters on which the external forces depend (e. g. 
the components of an external electric field). Forces and 
potentials will be supposed smooth, i.e. analytic, in their 
variables aside from possible impulsive elastic forces de- 
scribing shocks; the forces Fj will be supposed to vanish 
for 4> = 0. The impulsive forces are allowed here to 
model possible shocks with the walls of the container So 
or between hard core particles. 

Examples of deterministic reservoirs, [4], are forces ob- 
tained by imposing a nonholonomic constraint via some 
ad hoc principle, like the Gauss' principle, [5] (appendix 
9.A4), [6]. A different example will be extensively dis- 
cussed below. 

In general the forces i?i can be considered a set of 



deterministic "thermostat forces" if a further property 
holds: namely that the system evolves according to Eq. 
(1) towards a stationary state. This means that for all 
(X , X ), except possibly for a set of zero phase space vol- 
ume, any smooth function /(X , X ) evolves in time so 
that, denoting St(X°,X°) the configuration into which 
(X°, X°) evolves in time t according to Eq. (1), then the 
limit 

T lim ^ j T f(S t (X°,X°))dt = Jf(zMdz) (2) 

exists and is independent of (X° , X° ) . The phase space 
probability distribution fi(dz) (here z = (Y, Y) denotes 
the coordinates of position and velocity of a generic point 
in phase space) is then called the SRB distribution for 
the system. It has to be stressed that the condition that 
thermostat forces be "effective" enough to impede an in- 
definite build up of the energy of the system is a strong 
condition which we will assume in the models discussed 
in this note. It imposes on the interaction potentials and 
on the thermostats conditions that are not well under- 
stood, although they seem empirically verified with the 
simplest choices of molecular potentials, [5]. 

The maps St will have the group property S t ■ S# — 
St+t' and the SRB distribution /i will be invariant under 
time evolution. The SRB distribution is said to describe 
a stationary state of the system; it depends on the param- 
eters on which the forces acting on the system depend, 
e.g. |S| (volume), ^ (strength of the forcings) and on the 
model of thermostat forces. The collection of SRB distri- 
butions obtained by letting the parameters vary defines 
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a nonequilibrium ensemble. 

II. A MODEL OF EXTERNAL THERMOSTATS 

Let us now discuss in more detail an important class 
of thermostats in the context of model Eq. (1). Imag- 
ine that the No particles in the container So interacting 
via the potential V(X°) = J2 t<3 <P(4 ~ *§) + Ej M^) 
(where Vo models external conservative forces like ob- 
stacles, walls, gravity, . . .) and subject to the external 
forces Fj are also interacting with M other systems E , 
of N a particles of mass m a , in containers E a contiguous 
to So. It will be assumed that E a n E a / = for a ^ a', 
a, a' = 0, . . . ,M. 

The coordinates of the particles in the a-th system E a 
will be denoted x", j = 1, . . . , N a , and they will interact 

with each other via a potential V a (X a ) = Efj ^a( x f ~ 
x°). Furthermore there will be an interaction between 
the particles of each thermostat and those of the system 
via potentials W a (X°,X a ) = £f=i w a (x? - x»), 

a = l,...,M. 

Again, potentials will be assumed to be either hard 
core or non singular and Vo to be at least such that it 
forbids existence of obvious constants of motion. 

The temperature of each E will be defined by the total 
kinetic energy of its particles i.e., k B being Boltzmann's 

constant, by setting K a = J2f=i ! m a(*") 2 = \{dN a - 
l)fcs7o, where d is the spatial dimension; the particles of 
the a-th thermostat will be kept at constant temperature 
by further forces The latter are defined by imposing 
constancy of K a via Gauss' least effort principle. This 
means equations of motion like 

M 

mx° = -a x o(V(X°) + ^iy a (X°,X a )) +F 4 (X°,#) 

a=l 

m a x° = -d x a (K(X a ) + W a (X°, X a )) - 0? (3) 
and an application of Gauss' principle yields = 

(dM^-i^gT *j d ~ aa *j wnere L a is the work per unit 
time done by the particles in E on the particles of E a 
and V a is their potential energy, [6]. Note that in the 
first of Eq. (3), the forces — d x oW a play the role of the 
thermostat forces #j in Eq. (1). 

The work L a appearing in the definition of can 

be naturally interpreted as heat Q a ceded, per unit 
time, by the particles in Eo to the thermostat E a (be- 
cause the "temperature" of E a is constant). If X = 
(X°, X 1 , . . . , X M ), the entropy creation rate can be nat- 
urally defined as 

M A 

a—i 

hence exo can be called (in model Eq. (3)) the average 
entropy creation rate. Its time average will be assumed 



(To.+ = ( (T o)srb 7^ 0" Note that now (-)srb 1S ^ ne aver ~ 
age with respect to the stationary measure for the whole 
system E +thermostats. 

We shall see that, in ample generality, oo+ > 0: the 
definition of entropy creation is "reduced" , here, to an 
"equilibrium notion" because what is being defined is 
the entropy increase of the thermostats, which have to be 
considered in equilibrium. No attempt is made to define 
the entropy of the stationary state. Nor any attempt is 
made to define the temperature of the nonequilibrium 
system in E [T a is the temperature of E a , not of Eo). 

In fact the above model is a realization of a Carnot 
machine: the machine being the system in Eo on which 
external forces work leaving the system in the same state 
(a special "cycle") but achieving a transfer of heat be- 
tween the various thermostats (in agreement with the 
second law, see Eq. (4), only if a ,+ > 0). 

Another observable that is convenient to introduce is 
the phase space contraction rate <r(X,X), defined as mi- 
nus the divergence of the equations of motion: in the 
model described by Eq. (3), this is the sum of the deriva- 
tives of the r.h.s. with respect to X°, X a , mX°, m a X a 
and it turns out to be 

<7(X,X) = V %— ^ = y_%-n. (5) 

a kBTa a kBTa 

Therefore there is a simple and direct relation be- 
tween the phase space contraction and the entropy cre- 
ation rate, [6]: they just differ for the total deriva- 
tive u d = E a ktr ' w h° se ti me average (u} SRB van- 
ishes. This implies in particular that the time aver- 
ages cr + = (cr) SRB is the same as o- ,+ an d in particu- 
lar is nonnegative, consistently with the interpretation of 
phase space volume contraction rate. The cr + > implies 
that phase space contracts in average and therefore the 
SRB distribution will give probability 1 to a zero volume 
set. Therefore if cr+ > the system is said dissipative. 

The usefulness of introducing the definition of a is that 
a fluctuation theorem has been established for its large 
deviations functional in the context of Anosov systems 
theory. The chaotic hypothesis allows to establish a con- 
nection between the fluctuation theorem (valid for the 
phase space contraction rate in dissipative Anosov sys- 
tems) and a fluctuation relation for a and oo in model 
Eq. (3), as explained in next sections. 

III. CHAOTIC HYPOTHESIS 

In equilibrium statistical mechanics the ergodic hypoth- 
esis plays an important conceptual role as it implies that 
the motions have a SRB statistics and that the latter co- 
incides with the Liouvillc distribution on the energy sur- 
face. A role analogous to the ergodic hypothesis has been 
proposed for the chaotic hypothesis, [7] : which states that 

A chaotic system can be regarded as an Anosov system to 
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the purpose of computing the time averages of (smooth) 
observables. 

This means that the attracting set of a chaotic system, 
physically defined as a system with at least one posi- 
tive Lyapunov exponent, can be regarded as a smooth 
compact surface fi on which motion is highly unstable 
(uniformly hyperbolic) and transitive (there is a dense 
trajectory). For a mathematically precise definition of 
Anosov system we refer to [8]. 

We stress that the chaotic hypothesis concerns physi- 
cal systems: mathematically it is easy to find dynamical 
systems for which it does not hold. As it is easy (even 
easier) to find systems in which the ergodic hypothesis 
fails (e.g. harmonic lattices or black body radiation). 

Since physical systems are almost always not Anosov 
systems it is very likely that probing motions in extreme 
regimes will make visible the features that distinguish 
Anosov systems from non Anosov systems: much as it 
happens with the ergodic hypothesis. 

The ergodic hypothesis provides us with an expres- 
sion for the averages (as integrals over the normalized 
Liouvillc distribution on the energy surface) : likewise the 
Chaotic Hypothesis provides us with the existence and a 
formal expression for the averages (i.e. for the SRB dis- 
tribution), [8]. 

The interest of the hypothesis is to provide a frame- 
work in which properties like existence and formal expres- 
sion of an SRB distribution is a priori guaranteed, one 
can also say that the role of Anosov systems in chaotic 
dynamics is similar to the role of harmonic oscillators in 
the theory of regular motions. They are the paradigm 
of chaotic systems as the harmonic oscillators are the 
paradigm of order. Of course the hypothesis is only a be- 
ginning and one has to learn how to extract information 
from it, as it was the case with the use of the Liouville 
distribution once the ergodic hypothesis guaranteed that 
it was the appropriate distribution for the study of the 
statistics of motions in equilibrium situations, [9] . 



IV. FLUCTUATION THEOREM 

As mentioned above, an important observable in the 
theory of Anosov systems is the phase space contraction 
rate o~(x), defined as minus the divergence of the equa- 
tions of motion, computed in x G f2, where fi is the phase 
space of the system. A rather general result holds if the 
system is Anosov, dissipativc (cr + = (o~) SRB > 0), and 
furthermore reversible in the sense that there is an isom- 
ctry I of phase space such that IS t — S- t I for all t G R. 
Define the dimensionless phase space contraction 



p(x] 



T Jo 



o-(S t x) 



dt 



(6) 



then the SRB average of p is 1 (by definition) and there 
exists p* > 1 such that the probability P T of the event 



p G [a,b] with [a, b] C (— p*,p*) has the form 

P T (pE [a,b]) = const e ^max pe[a , 6] C(p)+0(l) ( 7 ) 

with C(p) analytic in (—p*,p*), [10, 11]. The function 
£(p) can be conveniently normalized to have value at 
p = 1 (i.e. at the average value of p). 

In Anosov systems which are reversible and dissipative 
a general symmetry property, called the fluctuation theo- 
rem and reflecting the reversibility symmetry, yields the 
parameterless relation, [7, 12], 



C(-p) = C(p) - P<t+ pz(-p*,p*) 



(8) 



where (—p*,p*), p* > 1, is the largest domain of defini- 
tion of (; it can be shown that £ is analytic on the whole 
(— p*,p*)- This relation is interesting because it has no 
free parameters. The relation was discovered in a simula- 
tion of a shear flow and it was suggested that it should be 
related to time reversal symmetry and to Ruelle's ideas 
on turbulence, [13]. A more informal (but imprecise) way 
of writing Eq. (7), (8) is 



Pr( P ) 

pA-p) 



_ e T PCT+ +0(l) 



for allpG (~p*,p*) (9) 



where P T (p) is the probability density of p. An interest- 
ing consequence of Eq. (9) is that (e _rpcr +) SRB = 1 in 
the sense that - log {e~ Tpa+ ) SRB T ^oo> 0- 

Occasionally systems with singularities have to be con- 
sidered: in such cases the relation Eq. (8) may change 
in the sense that the £(p) may be not analytic; one then 
expects that the relation holds in the largest analyticity 
interval symmetric around the origin. In various cases 
considered in the literature such interval appears to con- 
tain the interval (—1, 1) and sometimes this can be proved 
rigorously. For instance in simple, although admittedly 
special, examples of systems close to equilibrium, [9]. 

The Eq. (8), (9) is the first representative of con- 
sequences of the reversibility and chaoticity hypothe- 
ses. For instance given F\,...,F n arbitrary observ- 
ables which are (say) odd under time reversal / (i.e. 
F(Ix) = —F(x)) and given n functions t G [— J, J] — > 
ViWi 3 = one can ask which is the proba- 

bility that Fj(Stx) "closely follows" the pattern (fj(t) 
and at the same time - C (jl - SeX > d6 has value p. Then 

T JO CT+ r 

calling P T (F\ <~ <pi,...,F n ~ f n7 p) the probability 
of this event, which we write in the imprecise form 
corresponding to Eq. (9) for simplicity, and defining 

def 

Iifij(t) = — tpj(—t) it is 



Pr(Fl 



,F n 



Vn,p) 



P T (Fi ~ Ilfii, . . . , F n - Itp n , -p) 



(10) 



for p e (— p*,p*), which is remarkable because it is pa- 
rameterless and at the same time surprisingly indepen- 
dent of the choice of the observables Fj . The relation Eq. 
(10) has far reaching consequences, [5]. 
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Eq. (10) can be read as follows: the probability that 
the observables Fj follow given evolution patterns ipj con- 
ditioned to entropy creation rate po~+ is the same that 
they follow the time reversed patterns if conditioned to 
entropy creation rate — pcr+. In other words to change 
the sign of time it is just sufficient to reverse the sign of 
the average phase space contraction p (and we shall see 
that in our model this would amount at changing the sign 
of the entropy creation rate): no "extra effort" is needed. 

V. FLUCTUATION RELATION 

Given a chaotic system for which the Chaotic Hypoth- 
esis can be regarded as valid it is expected that the Fluc- 
tuation Relation (FR) holds: i.e. that one can define 
((p) and that the symmetry Eq. (8) holds. This is an 
important check that can be performed on the statistical 
properties of a stationary nonequilibrium state. 

It is however also important to know whether the quan- 
tity p and its fluctuations describe some interesting fea- 
ture of the dynamical system. The model Eq. (3), with 
a defined as in Eq. (5), provides an indication on the 
path to follow in the quest of an interpretation of the 
Fluctuation Relation [6, 14]. 

The remarkable property is that if we accept the 
Chaotic Hypothesis for the model Eq. (3) (i.e. for 
the whole system £ +thermostats) and we choose the 
parameters $ and T a in such a way that er + > 0, 
it is expected that FR holds for p: and this has an 
immediate physical interpretation because if we write 

Po(x) d = \ Jq a °j J g t { X ' > dt then, making use of the property 
that p — po is the variation of u in time r divided by the 
elapsed time r: 

u(S t x) - u(x) 
p(x) = po(x) + 

(70+ = {(?o)sRB = i a )sRB = a + ( n ) 

we find that in the limit t — > oo p and po have the same 
large deviation rate ((p). 

Of course the thermostats are "large" (we are even 
neglecting 0(N t f 1 )) and therefore the energies V a as well 
as u can be very large, physically of order 0(^2 a N a ). 
This means that the time t that has to pass to see the 
fluctuation relation Eq. (8) free of the 0(t _1 N a ) 
corrections can be enormous (possibly on astronomical 
scale in "realistic" cases). 

This is a situation similar to the one met when consid- 
ering systems with unbounded stochastic forces, [15], or 
with singular or nearly singular forces and thermostats 
of isokinetic type, [9]. In the present case we see that the 
really interesting quantity is the quantity p : which is 
the entropy creation rate and it is a boundary term un- 
affected by the size of the thermostats. Therefore if one 
considers and measures only p rather than p one not only 
performs a physically meaningful operation (i.e. measur- 
ing the average entropy creation rate) but also one can 



access the large deviation rate £(p) and check the Fluc- 
tuation Relation symmetry on a time totally unrelated 
to the thermostats size. Furthermore the more general 
relations like Eq. (10) can be naturally extended. 



VI. A MODEL FOR GRANULAR MATERIALS 

The current interest in granular materials properties 
and the consequent availability of experiments, e.g. [16], 
suggests trying to apply the above ideas to derive possible 
experimental tests. 

The main problem is that in granular materials colli- 
sions are intrinsically inelastic. In each collision particles 
heat up, and the heat is subsequently released trough 
thermal exhange with the walls of the container, sound 
emission (if the experiment is performed in air), radia- 
tion, and so on. If one still wants to deal with a reversible 
system, such as the one we discussed in the previous sec- 
tions, one should include all these sources of dissipation 
in the theoretical description. Clearly, this is a very hard 
task, and we will not pursue it here. 

A simplified description of the system consists in ne- 
glecting the internal degrees of freedom of the particles. 
In this case the inelastic collisions between point parti- 
cles will represent the only source of dissipation in the 
system. Still the chaotic hypothesis is expected to hold, 
but in this case the entropy production is strictly positive 
and there is no hope of observing a fluctuation relation, 
see e.g. [17], if one looks at the whole system. 

Nevertheless, in the presence of inelasticity, tempera- 
ture gradients are present in the system [16, 18, 19], so 
that heat is transported through different regions of the 
container. Then one can try to represent the processes of 
heat exchange between different regions of the system by 
the model we described above: and assuming that, un- 
der suitable conditions, the inelasticity of the collisions 
can be neglected, one can hope to observe a fluctuation 
relation for a (suitably defined) entropy production rate. 
This would be an interesting example of "ensemble equiv- 
alence" in nonequilibrium [5]: we will discuss this possi- 
bility in detail in the following. 

As a model for a granular material let S be a container 
consisting of two flat parallel vertical walls covered at 
the top and with a piston at the bottom that is kept 
oscillating by a motor so that its height is 

z(t) = Acosujt (12) 

The model can be simplified introducing a sawtooth mov- 
ing piston as in [19], however the results should not de- 
pend too much on the details of the time dependence of 
z(t). The container £ is partially filled with millimitcr 
size balls (a typical size of the faces of £ is 10 cm and the 
particle number is about 256): the vertical walls are so 
close that the balls almost touch both faces so the prob- 
lem is effectively two dimensional. The equations of mo- 
tion of the balls with coordinates (x{, Z{), i — 1,...,N, 
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Zi > z(t), are 

TTlXi — fx,i 

mzi = f z .i -mg + mS(zi- z(t))2(z(t) - Zi) (13) 

where ra=mass, g=gravity acceleration, and the colli- 
sions between the balls and the oscillating base of the 
container are assumed to be elastic [19] (eventually in- 
elasticity of the walls can be included into the model 
with negligible changes [17]); is the force describing 
the particle collisions and the particle- walls collisions. 

The force f, has a part describing the particle collisions: 
this is not necessarily elastic and in fact we will assume 
that the particle collisions are inelastic, with restitution 
coefficient a < 1. A simple model for inelastic collisions 
with inelasticity a (convenient for numerical implemen- 
tation) is a model in which collisions take place with the 
usual elastic collision rule but immediately after the ve- 
locities of the particles that have collided is scaled by a 
factor so that the kinetic energy of the pair is reduced by 
a factor 1 - a 2 [17-19]. 

We look at the stationary distribution of the balls: the 
simplest experimental situation that seems accessible to 
experiments and simulations is to draw ideal horizontal 
lines at heights hi > hi delimiting a strip So in the 
container and to look at the particles in So as a ther- 
mostatted system, the thermostats being the regions Si 
and S 2 at heights larger than h\ and smaller then hi, 
respectively. 

After a stationary state has been reached, the average 
kinetic energy of the particles will depend on the height 
z, and in particular will decrease on increasing z. Given 
the motion of the particles and a time interval t it will 
be possible to measure the quantity Q 2 of (kinetic) en- 
ergy that particles entering or exiting the region S from 
below (the "hotter side" ) carry out of So and the analo- 
gous quantity Q\ carried out by the particles that enter 
or exit from above (the "colder side"). 

If Tj, i = 1,2 are the average kinetic energies of the 
particles in small horizontal corridors above and below 
So, we see that there is a connection between the model 
of granular material, Eq. (13), and the model Eq. (3) 
discussed above. Still, model Eq. (13) cannot be exactly 
reduced to model Eq. (3) because of the internal dissipa- 
tion induced by the inelasticity a and of the fact that the 
number of particles in S depends on time, as particles 
come and go in the region. 

The reason for considering a model of granular material 
that is not in the class of models Eq. (3) is that Eq. (13) 
has a closer connection with the actual experiments [16] 
and with the related numerical simulations. Moreover, 
under suitable assumptions, that can be expected to hold 
on a specific time scale, the stationary state of Eq. (13) 
is effectively described in terms of the stationary state of 
Eq. (3), as discussed below. 

Note that real experiments cannot have an arbitrary 
duration, [16]: the particles movements are recorded by a 
digital camera and the number of photograms per second 
is of the order of a thousand so that the memory for the 



data is easily exhausted as each photogram has a size of 
about a 1Mb in current experiments. The same holds for 
numerical simulations where the accessible time scale is 
limited by the available computational resources. 

Hence each experiment lasts up to a few seconds start- 
ing after the system has been moving for a while so that a 
stationary state is reached. The result of the experiment 
is the reconstruction of the trajectory in phase space of 
each individual particle inside the observation frame, [16]. 

In order for the number of particles Nq in So to be ap- 
proximately constant for the duration of the experiment, 
the vertical size (hi — hi) of So should be chosen large 
compared to (Dt) 1 / 2 , where t is the duration of the ex- 
periment and D is the diffusion coefficient. Note that we 
are assuming that the motion of the particles is diffusive 
on the scale of So- In the low density case the motion 
could be not diffusive on the scale of So: then we would 
not be able to divide the degrees of freedom between the 
subsystem and the rest of the system and moreover the 
correlation length would be comparable with (or larger 
than) the size of the subsystem S . This would com- 
pletely change the nature of the problem and violations 
to FR could possibly be observed [20, 21]. 

Given the remarks above and if 

(1) we accept the chaotic hypothesis, 

(2) we assume that the result of the observations would 
be the same if the particles above So and below So were 
kept at constant total kinetic energy by reversible ther- 
mostats (e.g. Gaussian thermostats), [5, 22, 23], 

(3) we neglect the dissipation due to inelastic collisions 
between particles in So, 

(4) we neglect the fluctuations of the number of particles 
in S , 

(5) we suppose that there is dissipation in the sense that 

lim -(^ + ^)>0, (14) 

we expect the analysis of Sect. 5 to apply to model 
Eq. (13). 

Note that chaoticity is expected at least if dissipation is 
small and evidence for it is provided by the experiment in 
[16] which indicates that the system evolves to a chaotic 
stationary state in which dissipation occurs. Dissipation 
due to internal inelastic collisions will be negligible (to 
the purpose of checking a FR for cr ) only on a specific 
time scale, as discussed below. 

Accepting the assumptions above, we then predict that 
a Fluctuation Relation is satisfied, see Eq. (8) and (9), 
for fluctuations of 

P= T^{J\ + W (15) 

in the interval (— p*,p*) with p* equal (at least) to 1. 

The latter is therefore a property that seems accessi- 
ble to simulations as well as to experimental test. Note 
however that it is very likely that the hypotheses (2)- (4) 
above will not be strictly verified in real experiments, see 
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the discussion in next section, so we expect that the anal- 
ysis and interpretation of the experimental results will be 
non trivial. Nevertheless, the test would be rather strin- 
gent. 



VII. RELEVANT TIME SCALES 

The above analysis assumes the existence of (at least) 
two time scales. One is the "equilibrium time scale", 9 e , 
which is the time scale over which the system evolving 
at constant energy, equal to the average energy observed, 
would reach equilibrium in absence of friction and forc- 
ing. An experimental measure of e would be the decor- 
relation time of self-correlations in the stationary state, 
and we can assume that e is of the order of the mean 
collision time. Note that e also coincides with the time 
scale over which finite time corrections to FR become ir- 
relevant [24] : this means that in order to be able to mea- 
sure the large deviations functional for the normalized 
entropy production rate p in Eq. (15) one has to choose 
t >• 6 e , see also [25] for a detailed discussion of the first 
orders finite time corrections to the large deviation func- 
tional. A second time scale is the "inelasticity time scale" 
6d, which is the scale over which the system reaches a 
stationary state if the particles are prepared in a random 
configuration and the piston is switched on at time t = 0. 
Possibly a third time scale is present: the "diffusion time 
scale" Od which is the scale over which a particle diffuses 
over the size of So- The analysis above applies only if 
the time t in Eq. (15) verifies e <C t <C Od, Od (note 
however that the measurement should be started after a 
time ^ Od since the piston has been switched on in order 
to have a stationary state); in practice this means that 
the time for reaching the stationary state has to be quite 
long compared to e . In this case friction is negligible for 
the duration of the measurement if the latter is between 
e and min(0D,0d)- In the setting we consider, the role 
of friction is "just" that of producing the nonequilibrium 
stationary state itself and the corresponding gradient of 
temperature: this is reminiscent of the role played by 
friction in classical mechanics problems, where periodic 
orbits (the "stationary state") can be dynamically se- 
lected by adding a small friction term to the Hamilton 
equations. Note that, as discussed below, the tempera- 
ture gradient produced by friction will be rather small: 
however smallness of the gradient does not affect the "FR 
time scale" over which FR is observable [24]. 

If internal friction is not negligible (that is if t > Od) 
the problem would change nature: an explicit model 
(and theory) should be developed to describe the trans- 
port mechanisms (such as radiation, heat exchange be- 
tween the particles and the container, sound emission, 
...) associated with the dissipation of kinetic energy and 
new thermostats should be correspondingly introduced. 
The definition of entropy production should be changed, 
by taking into account the presence of such new ther- 
mostats. In this case, even changing the definition of en- 



tropy production it is not clear whether FR should be sat- 
isfied: in fact internal dissipation would break the time- 
reversibility assumption and, even accepting the Chaotic 
Hypothesis, nothing guarantees a priori the validity of 
FR. 

The validity of e <C t <C 0d,0o is not obvious in exper- 
iments. A rough estimate of Od can be given as follows: 
the phase space contraction in a single collision is given 
by 1 — a. Thus the average phase space contraction per 
particle and per unit time is a+.d = (1 — a)/0 e , where 
l/0 e is the frequency of the collisions for a given parti- 
cle. It seems natural to assume that Od is the time scale 
at which (J+^0d becomes of order 1: on this time scale 
inelasticity will become manifest. Thus, we obtain the 
following estimate: 

d ~ j^—Oe (16) 
1 — a 

In real materials a < .95, so that Od can be at most of 
the order of 200 e . Nevertheless it is possible that this 
is already enough to observe a Fluctuation Relation on 
intermediate times. 

The situation is completely different in numerical sim- 
ulations where we can play with our freedom in choosing 
the restitution coefficient a (it can be chosen very close 
to one [17-19], in order to have Od S> e ) and the size of 
the container So (it can be chosen large, in order to have 
D » e ). 

To check the consistency of our hypotheses, it has to be 
shown that it is possible to make a choice of parameters 
so that Od and Od are separated by a large time window. 
Such choices are possible, as discussed below. 

If (5 = h\ — hi is the width of So, e = 1 — a, 7 is the tem- 
perature gradient in So, and D the diffusion coefficient 
the following estimates hold: 

(a) e = 0(1) as it can be taken of the order of the in- 
verse collision frequency, which is O(l) if density is con- 
stant and the forcing on the system is tuned to keep the 
energy constant as e — > 0. 

(b) d = 0eO(e- v ) as implied by Eq. (16). 

(c) Od = O(jj) = 0(S 2 ) because D is a constant (if the 
temperature and the density are kept constant). 

(d) 7 = 0(^/e), as long as 5 <C e -1 / 2 . In fact if the 
density is high enough to allow us to consider the gran- 
ular material as a fluid, as in Eq. (5) of Ref.[19], the 
temperature profile should be given by the heat equation 
V 2 T + csT = with suitable constant c and suitable 
boundary conditions on the piston (T = T ) and on the 
top of the container (VT = 0). This equation is solved by 
a linear combination of const e ± ^/^ z ; which has gradients 
of order 0(y/e), as long as 5 <C 1/y/e and the boundaries 
of So arc further than 0(1/ ^/e) from the top. 

Now, if we choose <5 = e^ 13 , with [3 < \, and we 
take e small enough, we have e <C min{#<j,0£>} and 
S <C 0(e~2) ; as required by item (d). 
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Remark: The entropy creation rate due to heat transport 
into S , in presence of a temperature gradient 7, is 
given by er + = 0(7 2 <5) = 0(eS) because the temperature 
difference is 0(jS) and the energy flow through the 
surface is of order 0(7) (with 7 = 0(y/e), see item 
(d)). The order of magnitude of a + is not larger then 
the average amount ad of energy dissipated per unit 
time in S divided by the average kinetic energy T (the 
latter quantity is of order O(0~ 1 s5) because, at constant 
density, the number of particles in S is 0(6)); however 
the entropy creation due to the dissipative collisions in 
So has fluctuations of order 0(e5^ ) because the number 
of particles in S fluctuates by 0(5? ). This is consistent 
with neglecting the entropy creation inside the region 
S due to the inelasticity in spite of it being of the same 
order of the entropy creation due to the heat entering 
S from its upper and lower regions. 



This argument supports the proposal that in numeri- 
cal simulations it will be possible to test our ideas by a 
suitable choice of the parameters. We expect that other 
choices will be possible: for instance in the high-density 
limit it is clear that 6d S> 6 e because the diffusion coef- 
ficient will become very small. To what extent this can 
be applied to experiment remains an open question. 



VIII. REMARKS AND CONCLUSIONS 

(1) The model can be given further structure by adding 
a non conservative forcing acting on the particles in the 
region S : the same relations would follow (in partic- 
ular the Fluctuation Relation) if the forced equation of 
motion are still reversible; see [26] for a (stochastic) ex- 
ample. The above will not hold in general if the forcing 
is not reversible, e.g. if the inelasticity of the collisions 
inside S cannot be neglected, see below. 

(2) An explicit computation of the large deviation func- 
tion of the dissipated power, in the regime t ^> 8d [i.e. 
when the dissipation is mainly due to inelastic collisions) 
recently appeared in [27]. However in the model only 
the dissipation due to the collisions was taken into ac- 
count: so that it is not clear how the heat produced in 
the collisions is removed from the system, see the discus- 
sion above. It turned out that in this regime no negative 
values of p are observed so that the FR cannot hold. This 
is interesting and expected on the basis of the consider- 
ations above. It is not clear if, including the additional 
thermostats required to remove heat from the particles 
and prevent them to warm up indefinitely, the Fluctua- 
tion Relation is recovered. However this problem is dif- 
ferent from the one discussed here and we leave it for 
future investigation. 

(3) There has also been some debate on the interpreta- 
tion of the experimental results of [16]. In [17] a sim- 
plified model, very similar to the one discussed above, 
was proposed and showed to reproduce the experimental 



data of [16]. The prediction of the model is that the FR is 
not satisfied. Note however that the geometry considered 
in [16, 17] is different from the one considered here: the 
whole box is vibrated, so that the the temperature profile 
is symmetric, and a region So in the center of the box 
is considered. Heat exchange is due to "hot" particles 
extering S (Q+) and "cold" particles exiting S (Q-). 
One has Q = Q + + Q_ 7^ because of the dissipation in 
S and 



Q+ ^ Q- 



= 



(17) 



where T + is the temperature outside S and T_ is the 
temperature inside S , see [17]. Thus the dissipation due 
to heat exchange between the region S and the regions 
outside So vanishes and the only dissipation is due to 
the inelasticity of the collisions. In this regime, again, 
the FR is not expected to hold if the thermostat dissi- 
pating the heat produced in the collisions is not included 
in the model, see above: it is an interesting remark of 
[17] that partially motivated the present work. We be- 
lieve that different experiments can be designed in which 
the dissipation is mainly due to heat exchanges and the 
inelasticity is negligible, as the one we proposed above. 
The main difference between the experiment we proposed 
and the one in [16] is that the geometry of the box is such 
that Q1/T1 + Q2/T2 > 0: in this situation the dissipation 
due to inelastic collisions should be negligible as long as 

t<.e d . 

(4) Even in situations in which the dissipation is entirely 
due to irreversible inelastic collisions between particles, 
such as the ones considered in [17, 27], the Chaotic Hy- 
pothesis is expected to hold, and the stationary state to 
be described by a SRB distribution. In these cases the 
failure of the Fluctuation Relation is not in contradiction 
with the chaotic hypothesis, due to the irreversibility of 
the equations of motion. 

Conclusions We showed that in a paradigmatic class of 
mechanical models of thermostatted systems the phase 
space contraction is an interesting quantity. In large sys- 
tems in contact with thermostats it may consist of a sum 
of two quantities, the first with the interpretation of en- 
tropy creation rate and the second extremely large but 
equal to a total derivative. 

Its fluctuation properties, while asymptotically for 
large times determining (actually being identical to) the 
fluctuation properties of the entropy creation rate, hence 
implying the fluctuation relation, may require very long 
time to be freed of finite time corrections. But at the 
same time the study of the fluctuation properties of the 
physical quantity defined by the entropy creation rate can 
be used to determine the large deviations of the phase 
space contraction. The latter, having the same large de- 
viation rate, must obey the Fluctuation Relation which 
therefore becomes observable even if the system is in con- 
tact with large (or even infinite) thermostats. 

The analysis leads to propose concrete experimental 
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tests as well as tests based on simulations in the context 
of granular materials. The models naturally introduced 
for the description of granular materials experiments are 
not in the same class of models for which a relation be- 
tween phase space contraction and entropy production 
rate was previously discussed. The previous analysis can 
be applied to granular materials only under suitable as- 
sumptions, verified on a specific time scale. A fluctuation 
relation for an entropy production rate measured on the 
same time scale is predicted. 
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